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ABSTRACT 

We describe a model-independent method of assessing the uncertainties in 
cross-correlation lags determined from AGN light curves, and use this method to 
investigate the reality of lags between UV and optical continuum variations in 
well-studied AGNs. Our results confirm the existence of such lags in NGC 7469. We 
find that the continuum variations at 1825 A, 4845 A, and 6962 A follow those at 
1315 A by 0.22ig;ii days, 1.25±g || days, and 1.84±{];{g days, respectively, based on 
the centroids of the cross-correlation functions; the error intervals quoted correspond 
to 68% confidence levels, and each of these lags is greater than zero at no less than 
97% confidence. We do not find statistically significant interband continuum lags in 
NGC 5548, NGC 3783, or Fairall 9. Wavelength-dependent continuum lags may be 
marginally detected in the case of NGC 4151. However, on the basis of theoretical 
considerations, wavelength-dependent continuum lags in sources other than NGC 7469 
are not expected to have been detectable in previous experiments. We also confirm the 
existence of a statistically significant lag between X-ray and UV continuum variations 
in the blazar PKS 2155-304. 



Subject headings: galaxies: active — galaxies: Seyfert — methods: data analysis 



1. Variability in Active Galactic Nuclei 

Over the last ten years, a number of intensive monitoring experiments on active galactic 
nuclei (AGNs) have been carried out (for reviews, see Netzer & Peterson 1997 and Peterson 
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1993). While the primary purpose of many spectroscopic monitoring campaigns has been to 
determine the response of the broad emission lines to continuum variations and thus determine the 
structure and kinematics of the line-emitting gas through the process of "reverberation mapping" 
(Blandford & McKee 1982), these campaigns have also provided an opportunity to search for time 
delays between different continuum bands. Such time delays are expected to exist if the continuum 
flux in one waveband is reprocessed continuum emission from another waveband, either directly 
through irradiation or indirectly through viscous processes of the emitting plasma. The standard 
accretion-disk model for AGNs may exhibit both. 

Recently, Wanders et al. (1997) found evidence for such time delays in the UV spectra of the 
Seyfert 1 galaxy NGC 7469 obtained with the International Ultraviolet Explorer (WE) spacecraft. 
It was found that the continuum flux variations around 1700-1800 A lag behind those at 1315 A by 
0.3 ± 0.07 d. This result was subsequently supported by contemporaneous optical spectra (Collier 
et al. 1998) that showed that the continuum variations around 4825 A lag behind the 1315 A 
variations by 1.2 ± 0.3 d, and those around 6925 A lag by 1.7 ± 0.7 d. The increase in time lag with 
wavelength in NGC 7469 is strong evidence for continuum reprocessing models of AGNs in which 
the longer-wavelength photons are reprocessed shorter-wavelength photons originating closer to 
the central source. 

A legitimate question to ask at this point is why such an effect has not been previously 
reported in the literature, even though other AGNs have been monitored in a similar fashion? 
This question is of particular interest in the case of the well-studied galaxy NGC 4151, which 
was monitored at an even higher sampling rate than NGC 7469 during a 10-day campaign in 
1993 (Crenshaw et al. 1996; Kaspi et al. 1996; Warwick et al. 1996; Edelson et al. 1996). Other 
Seyfert galaxies that have been well-monitored simultaneously in the UV and optical are NGC 
5548 (Clavel et al. 1991; Peterson et al. 1991, 1992; Korista et al. 1995), NGC 3783 (Reichert 
et al. 1994; Stirpe et al. 1994), and Fairall9 (Rodriguez-Pascual et al. 1997; Santos-Lleo et al. 
1997). None of these studies found a significant lag between different continuum waveband flux 
variations. Why were such lags found for NGC 7469, but not for the others? 

The answer to this question is not trivial. The detection of such lags is dependent on 
the sampling characteristics of the light curves, as well as on the variations of the light curves 
themselves, i.e., the auto-correlation function (ACF) of the light curve. In any event, the formal 
statistical significance of a time-lag detection through a cross-correlation analysis is hard to assess. 
There is no generally agreed upon way to estimate the errors in cross-correlation lags and thus 
attach a level of significance to a time-lag detection. Indeed, interband continuum lags have been 
reported in previous campaigns, but the detections were not thought to be statistically significant. 

In this paper, we will introduce what we believe is a conservative model-independent approach 
to estimating the uncertainties in cross-correlation lag determinations (§ 2) and use it to re-examine 
existing data sets (§ 3). We will then discuss the implications of these results (§ 4) and present our 
conclusions (§5). 
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2. Estimation of Cross-Correlation Uncertainties 

In studies of AGN time series, we are faced with the following general situation: we have two 
light curves with fluxes {A±,A2, . . . , An} and {B\, B2, . . . , Bm}, that are often irregularly sampled 
at times {ti, £2, • • • , i?v} and {t[, t' 2 , . . . ,t' M }, respectively. The two light curves have ./V and M 
observed fluxes Fj with corresponding estimates of the la measurement errors {01,02, . .. ,ajy} 
and {/3i,/?2, • • • ,Pm}- We then cross-correlate the two light curves. The quantity we seek is 
the lag (whether it be the lag of the peak or the centroid of the CCF is irrelevant at this 
point), and its uncertainty. In general, we recommend using the centroid of the cross-correlation 
function T cen t rather than the location of the peak T pea k- Koratkar & Gaskell (1991) and Penston 
(1991) have shown that T cen t is a measure of the luminosity-weighted radius of the reprocessing 
region. However, as Robinson & Perez (1990) have shown, T pea k cannot be interpreted as a 
physical quantity in any simple manner. Besides being model-dependent, it also depends on the 
auto-correlation function of the light curves, as well as the sampling of the light curves. The 
results of their simulations showed that the peak lag was most sensitive to the inner regions of a 
reprocessing region (also as noted by Gaskell Sz Sparke 1986), and thus likely underestimates the 
actual size of the reprocessing region. 

2.1. Model-Dependent Monte Carlo Methods 

While there is no general agreement on how to estimate the uncertainties in cross-correlation 
lags, probably the method that has the widest acceptance is Monte Carlo (MC) simulations. In 
studies of AGN variability, these are generally carried out in the following fashion: first, one 
adopts a model light curve that is supposed to "drive" the other variations. Selection of a driving 
light curve is usually based on the shortest UV continuum wavelength observed, and it is assumed 
that the ionizing continuum (which is supposed to drive the emission-line variations) behaves in a 
similar fashion. The driving light curve is then convolved with a "transfer function" that produces 
a "responding" light curve, which most often represents the response of a broad emission line to 
the continuum variations. The driving and responding light curves are then "sampled" in a fashion 
that somehow mimics the real observations, and then the effects of observational uncertainties, 
both random and systematic, are included. These artificial data sets are then cross-correlated as 
if they were real data, and the cross-correlation lag (based either on the peak or centroid of the 
cross-correlation function) is recorded. This process is repeated a large number of times to build 
up a cross-correlation peak distribution (CCPD; see Maoz & Netzer 1989) which we emphasize 
is not necessarily a normal distribution. From the CCPD, the likelihood that a given lag falls in 
some particular range can be determined. 

The principal problem with this particular method is that the results are highly model- 
dependent — the uncertainty we derive is only valid to the extent that our models of the driving 
light curve and the transfer function are reasonable approximations to the true situation. 
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The behavior of the continuum is one of the principal sources of uncertainty in MC 
simulations. AGNs undergo irregular continuum variations whose origin is not known, but thought 
by some to have some relationship to accretion-disk instabilities or changes in the accretion rate 
(e.g., Mineshige & Shields 1990; Siemiginowska, Czerny, & Kostyunin 1996). AGN-like continuum 
variations can be parameterized by a power-density spectrum (PDS) of power-law form P(f) oc f a 
(where / is the temporal frequency in Hz or days -1 ) and index a w — 2 ± 0.5. A model continuum 
that bears a close resemblance to observed continua can therefore be generated from a PDS of 
the correct form by randomizing the phases of the Fourier components, performing the inverse 
Fourier transform, and then normalizing the continuum to yield the desired fractional variability 
(F var , as defined by Rodriguez-Pascual et al. 1997). For a limited data train, the variations in 
such model continua can show a very wide variety of structures. Indeed, for well-sampled high 
signal-to-noise ratio light curves, the principal source of variance among various MC realizations 
is often determined by differences in continuum light curve features, even when all the model 
continua are characterized by the same PDS power-law index and fractional variability. In such 
cases, MC simulations can in fact overestimate the uncertainty in cross-correlation lags. Certainly, 
however, such simulations are extremely valuable and valid in experiment design, since the 
particulars of the continuum behavior cannot be predicted. 

Determination of a model responding light curve is similarly fraught with uncertainty, not 
only because of the unknown geometry of the reprocessing region, but because the reprocessing 
physics needs to be modeled as well (regardless of whether we are considering reprocessing into 
emission lines or continuum radiation). 

We wish to consider, therefore, methods by which the uncertainties in cross-correlation lags 
can be estimated in a less model-dependent fashion that is simple to implement. 

2.2. Model-Independent Monte Carlo Methods 

Given a sufficiently long, continuously sampled light curve with error-free observations, it 
should be possible to determine a cross-correlation lag to arbitrary precision. The experimentally 
determined lag is not necessarily highly repeatable in time (even if the transfer function does not 
change with time) simply because the specific pattern of continuum variability can change with 
time; since the responding light curve is a convolution of the transfer function and the continuum 
light curve, the line response can also show pronounced differences in repeat experiments. It 
is simple to show (e.g., Penston 1991) that the cross-correlation function itself is a convolution 
of the transfer function and the driving continuum auto-correlation function, so the centroid or 
peak of the cross-correlation function can vary on account of differences in the driving continuum 
autocorrelation function. We therefore consider only sources of uncertainty that are introduced by 
deficiencies in the experimental data, and not the broader problem of attempting to determine 
limitations on the transfer function imposed by the measurements. 
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Two of the principal sources of uncertainty in cross-correlation lags are (a) flux uncertainties 
in individual measurements and (b) uncertainties associated with the observational sampling of 
the light curves, i.e., the intervals between observations and the duration of the experiment. We 
consider methods for assessing the uncertainties associated with each of these. 

In many AGN monitoring programs, uncertainties in measured fluxes can be a significant 
source of error. The importance of flux errors can be easily assessed through MC simulations. 
One just takes each real measurement Ai and alters it by adding a random noise contribution. We 
assume that the errors in fluxes are normally distributed, and we thus modify each flux by random 
Gaussian deviates based on the quoted error ccj for each datum. The modification of each data 
point is statistically independent from each of the others. In a single MC realization, each data 
point is modified, the CCF is computed, and the lag is recorded. Multiple realizations build up 
the CCPD. Comparing a large number of independent realizations should reveal that the average 
value of an individual point remains Ai, and the standard deviation should be c^. We will refer 
to this process as "flux randomization" (FR). This procedure is commonly used as part of MC 
simulations, as described above. 

Cross-correlation results can be highly sensitive to individual data points or even combinations 
of data points. We can, then, to some extent test the sensitivity of a cross-correlation result 
by considering only subsets of the original parent data set. This process is similar, but not 
identical, to a commonly used and powerful statistical technique known as "bootstrapping". 
The bootstrap method can be used to evaluate the significance of correlations based on limited 
data (see Press et al. 1992 for a brief description). For example, consider a set of observations 
{(xi, y±), (x2, 2/2), (%n, Vn)}- For a single realization, we draw from this sample N pairs 
of randomly selected points (xj,t/j) without regard to whether or not they have been previously 
selected. This results in exclusion of some data pairs, and counts others multiply. The correlation 
analysis is again performed, and the correlation coefficient is recorded. Multiple realizations lead 
to a mean and standard deviation for the correlation coefficient. 

Extension of this process to time series, in which one set of values is highly ordered, is not 
obvious. The temporal order of the data points must be preserved. But we can still select a 
subset of the original data points without regard to previous selection; how to incorporate these 
multiple data points into the analysis becomes somewhat problematic: one can either ignore the 
redundant data points (which effectively reduces the size of the selected sample, typically by 
a factor ~ 1/e ~ 0.37, which is the Poisson probability of not selecting any particular point), 
or one can increase the weight of the multiply selected data points. The latter can be used in 
cross-correlation methods that make use of the flux uncertainties, such as the discrete correlation 
function (DCF) method (Edelson & Krolik 1988) or the z-transformed discrete correlation function 
(ZDCF) method (Alexander 1997). Here we will use only the interpolated cross-correlation 
function (ICCF) method, which does not make use of the flux errors in the data, and therefore we 
will simply exclude redundant selections from the simulations. Thus, each realization is based on 
a randomly chosen subset of the original data points, and we will refer to this as "random subset 
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selection" (RSS). We argue that this method gives a fairly conservative estimate of the errors due 
to sampling: the effects of individual data points are accounted for statistically by random removal 
of such points, and each MC realization is based on a number of data points that is typically 37% 
smaller than the real data set, so the individual MC realizations are likely to produce a broader 
range of values than similar sets the same size as the original set. 

The two elements that we describe above are complementary in that they test the sensitivity 
of the cross-correlation results to flux errors (FR) and to sampling characteristics (RSS), and 
the two processes can be combined in a single simulation. This combination does not account 
for all possible sources of error, however, as we still cannot test what might have occurred in 
the gaps between observations. In all realizations, the data are restricted to some subset of 
the actual observation times {t±, • • • , ijv} and t' 2 , ■ ■ ■ , t' M }. Dependence on the times of 
observations can be tested only if one has a model for the flux behavior at other times — this 
leads to model-dependent simulations of the sort described above. Our proposed method should 
be valid, however, to the extent that the ICCF method is itself valid, i.e., the gaps between 
observations must be sufficiently small that low-order interpolation between them is a reasonable 
approximation to the true behavior of the light curve. 

We also caution the reader that the method as implemented here does not take into account 
systematic errors that may be operative in some cases. An obvious example is the problem 
of "correlated errors", which can arise when multiple measurements are made from a single 
spectrum, which is the case for many of the data discussed here. Suppose, for example, that 
the flux-calibration of a set of spectra has a random-error component. Then a time series of 
measurements of a even a constant source will show (presumably small) variations, with all 
measured spectral features varying similarly and in phase. Cross-correlation of these time series 
will yield a peak at zero lag. In the case of AGN variability studies, correlated errors have been 
found to introduce significant systematic effects only in one case, Akn 120 (Gaskell & Peterson 
1987). This was a result of the relatively small equivalent width of the [Oiii]A5007 emission 
line (which is commonly used for flux calibration in optical variability studies of AGNs) and the 
relatively large measurement uncertainties (typically ~8%) that could be obtained with pre-CCD 
era detectors. On the basis of simulations, Gaskell & Peterson (1987) concluded that the effects of 
correlated errors would be obviated by a modest improvement in the quality of the data, and this 
has indeed been borne out (Peterson et al. 1998). As a cautionary measure, we have examined 
some of the data sets used here for correlated errors. We have looked for correlations between the 
differences of pairs of measurements closely spaced in time, i.e., correlations between A(tj) — A(U) 
and B(tj) — B(t{), where A and B are measured from the same spectra. We find no evidence 
that correlated errors affect the data discussed here. In any event, we note that if correlated 
errors were in fact present, the artificial correlation at zero lag would tend to drive our measured 
cross-correlation lags to smaller values, thus masking the very effect we are attempting to measure. 
If correlated errors affected these data, our results would thus be conservative estimates of the 
probability that the interband lags are greater than zero. 
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2.3. Test Results 

In order to test our assertion that the combined FR/RSS method will yield conservative error 
estimates of cross-correlation lags, we will undertake a number of numerical tests based on simple 
models, which we construct as follows: first a driving continuum is constructed from a power 
spectrum of a given PDS index a and randomized phases. This process is described more fully 
by White & Peterson (1994). The spectrum produced is then normalized to a root-mean-square 
(rms) variability level F var that is consistent with experimental results. The driving continuum 
is then convolved with a specified transfer function to produce a responding light curve, which 
can then be renormalized to a value of F var that is consistent with what is observed for emission 
lines. These two light curves constitute "parent" light curves from which individual experimental 
realizations can be drawn. 

For each realization, we select according to some specified procedure TV" data points spanning 
some time span tN — h from both the driving and responding (which we will call respectively the 
"continuum" and " emission-line") light curves. In these tests, the continuum and emission-line 
fluxes are always paired, i.e., both time series have N points sampled at the same time; this 
is not a necessary restriction, but it mimics what occurs in nearly all real data sets. Once a 
sample of data points has been selected, the effects of random errors are included. Fractional 
flux deviations based on random Gaussian deviates with variances again chosen to be consistent 
with real data sets are then used to redistribute the sampled fluxes and thus simulate the effect 
of random flux errors. The N selected data points from each series are then cross correlated as if 
they were real data. A cross-correlation realization is deemed to have succeeded if the maximum 
correlation coefficient r max is large enough to conclude that the correlation between the two series 
is significant at a level of confidence greater than 95%. The time delay T pea k at which r max occurs 
and the centroid of the cross-correlation function T cen t, here based on all points with correlation 
coefficients in excess of 0-8r max , are recorded and used to build up distribution functions (CCPDs) 
for these parameters. After a large number of realizations (usually 500 to 2000), this distribution 
can be integrated to determine the probability that a given realization will yield a result above or 
below some specified value. We reiterate the important point made by Maoz & Netzer (1989) that 
these distributions are almost always non-normal and the standard deviation of the distribution is 
not a good characterization of the uncertainty in the lag. Unless otherwise noted, the uncertainties 
we will quote, computed directly from the CCPDs, will be ±At68, defined such that 68.27% of 
the realizations yield results between T me dian — At"68 and T me dian + At68> or equivalently, 15.87% 
of the realizations give values below T me dian — Atqs and an equal number above T me dian + Atqs; 
these limits correspond to la errors for a normal distribution. A sample CCPD for one of the tests 
performed below is shown in Fig. 1. 

In order to compare the results of our method with the results of model-dependent MC 
simulations, we select at random one of the realizations for a given model, and we independently 
attempt to estimate the uncertainties in the lag by FR/RSS, and compare these to the results of 
the model-dependent simulations. The results of a single sample test are shown in Table 1. For 
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the test described here (based on one of the White & Peterson 1994 models with a AT = 1 day 
interpolation grid and where the responding emission-line region is modeled as a thin spherical 
shell), the continuum PDS has index a = —2.5, and F var = 0.30. The BLR is characterized by 
R = 20 light days and asymmetry parameter A = (i.e., completely isotropic line emission). The 
assigned fractional error levels are 0.035 for the continuum and 0.037 for the line. Each sample 
consists of N = 40 data points, chosen at random within a 200-day span. Column (1) gives 
the parameter recorded, either the centroid T cen t or peak T pea k of the cross-correlation function. 
Column (2) gives the median value and uncertainties (i.e., ±At68) for the model-dependent MC 
simulation. Column (3) gives the cross-correlation lag {r cen t or T pea k) based on the single randomly 
selected realization and the FR/RSS-based uncertainty computed from that realization. For 
comparison, columns (4) and (5) give respectively the uncertainties based on RSS and FS alone. 
The results of this randomly chosen test are representative in several respects. First, as expected, 
the error estimates based on either RSS or FR alone are generally smaller than the those for the 
combined FR/RSS. Second, the expected value of the lag from the model (20 days) is within the 
error limits for each of the methods. Third, the combined FR/RSS uncertainties are somewhat 
larger than the uncertainties based on the model-dependent MC simulations, i.e., they seem to 
provide a conservative estimate of the uncertainties. However, the uncertainties based on either 
FR alone or RSS alone are usually smaller than the MC-based uncertainties, and we will therefore 
not consider these further. 

As we mentioned earlier, specification of the continuum with only the PDS index a and the 
fractional variability F var can produce a wide variety of light curves. In Table 2, we show the 
results from 10 nominally identical simulations, quoting both the model-dependent Monte Carlo 
results and the FR/RSS results based in each case on one randomly selected MD realization (Test 
1 is from Table 1). The entries in Table 2 show that T cen t and T pea k vary somewhat, depending on 
the precise continuum behavior that is sampled. Moreover, the uncertainties in these parameters 
can also vary. The FR/RSS uncertainty is always greater than or approximately equal to the 
model-dependent Monte Carlo uncertainty, and the true value of 20 days is virtually always within 
the quoted uncertainty range. 

In Table 3, we show the results for a similar series of tests on a thick-shell model with good 
temporal sampling (N = 50 observations randomly spaced on a grid of spacing AT = 0.5 days 
during a 100-day interval). The continuum used is characterized by a power-law PDS with a = 2 
and F var = 0.14. The specific thick-shell model chosen is the Wanders et al. (1995) for NGC 
5548, with an inner radius of 1 light day, and outer radius of 12 light days, and an emission-line 
asymmetry A = 0.6 (although we do not here consider the effects of an anisotropic continuum 
source, as did Wanders et al.). Both lines and continuum are assumed to be measured to 2.5%. 
The results are shown in Table 3. These tests show that lag can vary significantly depending on 
the details of the light curve, which are different for each of the 10 tests shown here. Again it 
is seen that the FR/RSS errors are similar to or slightly greater than the model-dependent MC 
errors, even though the errors themselves also vary from test to test. 
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In Table 4, we show more extensive comparisons of the model-dependent MC-based and 
FR/RSS uncertainties. The models shown in Table 4 are a subset of the thin-shell models that 
cover most of the parameter space considered by White & Peterson (1994). Column (1) gives the 
PDS index for the continuum model and column (2) gives the broad-line region radius (the transfer 
function for a thin spherical shell of radius R is a rectangular function that is non-zero only in 
the range to 2R/c, thus with centroid R/c). Column (3) gives the number of observational data 
points sampled from a 200-day span of data — the mean sampling intervals are thus about 5 
days for N = 40, 10 days for N = 20, and 20 days for N = 10. The cases N = 20 and N = 10 
constitute rather poor sampling by AGN monitoring standards, particularly in the case of the 
flatter PDS (a = —1.8) which has more rapid continuum variations. Columns (4) and (5) give 
the median value of the centroid and associated uncertainties for the complete model-dependent 
Monte-Carlo simulations (MC) and for the model-independent (FR/RSS) simulations, respectively. 
The corresponding values of T pea k are given in columns (6) and (7). These results show that 
in general the FR/RSS uncertainties are quite comparable or slightly more conservative than 
the model-dependent MC errors. As the sampling rate is decreased (i.e., as N decreases), the 
uncertainties in the cross-correlation results become very large, as can be seen in the Table. In the 
case of very poor sampling, the FR/RSS uncertainties are much larger than the MC errors because 
the typical FR/RRS realization is working with about 37% fewer points from an already poorly 
sampled light curve. Thus, in the limit of very poor sampling, the FR/RSS method fails in a 
desirable way by giving error estimates that are too large (in contrast, say, to the Gaskell-Peterson 
[Gaskell & Peterson 1987] analytic formula, as shown by White &; Peterson 1994). 

In summary, we conclude that FR/RSS gives reliable, conservative uncertainties for a wide 
variety of cross-correlation simulations. We now proceed to apply this method to estimate errors 
for measured interband continuum lags for a number of AGNs that have been simultaneously 
monitored in the UV and optical. 

3. Wavelength-Dependent Lags in Well-Studied Sources 

In Table 5, we show cross-correlation results for NGC 7469, where the UV data are from 
Wanders et al. (1997) and the optical data are from Collier et al. (1998). In addition to the cross- 
correlation lags T cen t and r pea k, we show estimates of the uncertainties in each of these parameters, 
based on FR/RSS Monte Carlo simulations as described in the previous section. We also show 
the probability, again based on the FR/RSS simulations, that the true lag is less than or equal to 
zero. As noted earlier, the FR/RSS uncertainties are conservative compared to model-dependent 
Monte-Carlo simulations, and thus the uncertainties reported here are somewhat larger than 
quoted in the original work (for example, for the 1315 A/1825 A cross-correlation, Wanders et 
al. find T cen t = 0.22 ± 0.07 days for the TOMSIPS-derived data, and we assign uncertainties of 
+0.12/ — 0.13 days). Even using the FR/RSS errors, the wavelength-dependent continuum lags 
are statistically significant (although not r pea k for the 1315 A/1825 A cross-correlation), though 
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not highly significant. 

A number of additional multiwavelength AGN monitoring campaigns have been carried out 
over the past years, and these should provide us with data suitable for searching for similar 
wavelength-dependent continuum lags in other AGNs. Several non-blazar AGNs have been 
monitored simultaneously in the optical and ultraviolet 1 , but with the exception of NGC 7469, 
these have yielded only upper limits on possible lags between continuum variations in the 
short-wavelength ultraviolet (usually at about 1300 A) and the optical (usually at about 5000 A). 
However, the upper limits quoted depend on the uncertainties in the lag determination, and these 
have not been estimated in a uniform fashion. Here we employ the FR/RSS methodology to 
re-evaluate the uncertainties in cross-correlation lags for these well-studied AGNs in a uniform 
fashion in order to investigate whether or not the NGC 7469 results may be atypical. We will 
also include for completeness results on one well-studied blazar. The AGNs we consider are the 
following: 

1. NGC 5548: We re-examine the results of two UV/optical campaigns, the first from 1989 
(Clavel et al. 1991, Peterson et al. 1991, 1992) and a more intensive campaign in 1993 
(Korista et al. 1995). In the earlier campaign, the IUE data covered an eight-month period 
from 1988 December to 1989 August, yielding measurements at 60 epochs with an average 
interval of about 4 days between observations. Concurrent ground-based spectra yielded 
optical continuum measurements at an average spacing of about 2.5 days. The 1993 program 
was anchored by Hubble Space Telescope (HST) FOS observations which were obtained 
once per day for 39 days. The program was preceded by a series of IUE observations 
(one every other day) which together yield a light curve consisting of 53 data points over 

a 76-day period. The concurrent ground-based program produced 105 optical continuum 
measurements during the same period. 

2. NGC 3783: This galaxy was monitored with IUE from 1991 December to 1992 July 
(Reichert et al. 1994), yielding 69 observations with an average spacing of 3.3 days. 
Concurrent ground-based observations (Stirpe et al. 1994) produced an optical light curve 
that was comprised of 72 observations during the same period with approximately the same 
average spacing as the IUE data. 

3. Fairall 9: IUE observations (Rodriguez-Pascual et al. 1997) were obtained between 1994 
April and 1994 December. Optical observations were obtained from the beginning of this 
program through the end of 1995 January (Santos-Lleo et al. 1997); however, the sampling 
became markedly worse late in the campaign, so we restrict our attention here to the data 
obtained through early October. The data considered here consist of 39 IUE measurements, 



1 All of the non-blazar AGNs studied here were monitored by the International AGN Watch consortium. The light 
curves used in these calculations can be obtained through the International AGN Watch site on the World-Wide Web 
at URL http://www.astronomy.ohio-state.edu/~agnwatch/. 
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with an average spacing of 4.2 days between epochs, and 25 optical measurements with an 
average spacing of slightly more than 6 days. 

4. NGC 4151: Prior to the NGC 7469 campaign, NGC 4151 represented the state-of-the-art 
in multiwavelength studies of AGN variability, with intensive observations obtained over 

a 10-day period in 1993. The UV data, with a mean spacing of about 0.05 day, are from 
Crenshaw et al. (1996) and the optical data, with a mean spacing of about 0.6 days, are from 
Kaspi et al. (1996). 

5. PKS 2155—304: This is an especially well-studied blazar that was monitored intensively 
in soft X-rays (Brinkmann et al. 1994), the ultraviolet (Urry et al. 1993), and the optical 
(Courvoisier et al. 1995) in 1991 November. Analysis of the light curves by Edelson et al. 
(1995) indicates that the variations in the soft X-rays precede those in the UV/optical by 
~2-3 hours. 

For each of these AGNs, we have cross-correlated the data described above using the 
interpolation cross-correlation as outlined by White & Peterson (1994), using only the parts 
of the light curves that overlap in time, and have estimated the uncertainties by using the 
FR/RSS method. The results for three of the galaxies are given in Table 6. In both NGC 5548 
monitoring campaigns, it was found that the optical continuum (at 5100 A) lagged behind the 
UV continuum (at about 1350 A) by T cen t ~ 1 day. The uncertainties computed by the FR/RSS 
method indicate that the detected lag is not statistically significant, as was previously concluded 
from model-dependent Monte Carlo simulations. Similar conclusions are reached for both NGC 
3783 and Fairall 9, which were less favorable for detection of such small lags, and thus have larger 
associated uncertainties. 

The case of NGC 4151 is somewhat more complicated, as shown in Table 7. The 
1275 A/1820 A, 1275 A/2688 A, and 1275 A/5125 A cross-correlations have already been computed 
by Edelson et al. (1996), including maximum likelihood error estimates Atml on r pea k', these 
values are given in column (3) of Table 7. It can be seen that the FR/RSS estimates are somewhat 
smaller for the 1275 A/1820 A and 1275 A/2688 A cross-correlations, but similar for 1275 A/5125 A. 
As concluded by Edelson et al., the lags for all of these bands relative to the 1275 A continuum are 
not statistically significant. However, the 1275 A/6925 A cross-correlation, which was not reported 
by Edelson et al. and is reported here for the first time, shows that the 6925 A continuum lags 
behind the 1275 A continuum by about 2.5 days, at almost the 90% confidence level. 

The UV/optical CCFs for NGC 4151 have a broad and noisy structure, and are not very 
well-defined (see Fig. 5 of Edelson et al. 1996). A significant lag between the UV and optical 
continuum is thus not obtained, with the possible exception of the 6925 A light curve. The large 
uncertainties are mainly a result of the very few optical observations and the fact that the optical 
continuum does not trace the highest UV variability frequencies (Edelson et al. 1996). 

Given that we have apparently detected a marginally significant lag relative to 1275 A for 
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only one waveband, it is worthwhile re-examining all wavelengths in this range, as the existence of 
a wavelength-dependent trend in the data may be more convincing, as indeed it was in the case of 
the UV data alone in NGC 7469 (see Fig. 6 of Wanders et al. 1997). We have done this for both 
the UV and optical data by making new measurements from the original spectra. 

All SWP UV spectra of NGC 4151 were rebinned into 20 A-wide bins, and all LWP UV were 
rebinned into 50 A-wide bins. The optical spectra from Wise Observatory were similarly averaged 
into 100 A-wide bins. For each bin, we construct a light curve by measuring all of the spectra in 
the time series, and we then cross-correlate these light curves with the 1275 A light curve (for the 
UV light data) curves and with the 5125 A light curve (for the optical light curves). We then find 
the centroid lags as a function of wavelength. The results are plotted in Fig. 2 for the UV data 
and in Fig. 3 for the optical data. For the UV data, Fig. 2 shows that the lags are consistent with 
a constant value of zero, independent of wavelength. For the optical data, this is not so obvious. 
Clearly, there appears to be an increase in lag with wavelength, at least for wavelengths longer 
than about 5100 A. For shorter wavelengths, the light curves are contaminated with emission-line 
features, which have a longer lag than the continuum (Kaspi et al. 1996), and these wavelengths 
are not representative of the true continuum. For the longer wavelengths, we see contamination 
by the Ha emission line, but other than that, there is very little line contamination. 

Although the origin of the continuum emission in blazars is almost certainly different than in 
other AGNs, we include the results on the blazar PKS 2155—304 for completeness and as another 
illustration of the efficacy of the FR/RSS method. The results are given in Table 8; note that the 
uncertainties in this Table are 90% confidence limits rather than 68% confidence limits for direct 
comparison with Table 2 of Edelson et al. (1995). Comparison of the FR/RSS results with those 
from several other methods reveals that the FR/RSS errors are only slightly larger than those 
based on the \ 2 minimization method, which Edelson et al. regard as the most powerful technique 
for these particular data. 

In summary, we find that in the cases of NGC 5548, NGC 3783, and Fairall 9, statistically 
significant wavelength-dependent continuum lags are not found, simply because the uncertainties 
in the lag determinations are relatively large, greater than about 1 day in each case. Only in 
the case of NGC 4151, which is somewhat less luminous than NGC 7469 and therefore might be 
expected to show more rapid variability, is the temporal sampling good enough that one might 
hope to see wavelength-dependent continuum lags. And indeed, there is at least weak evidence 
that this effect also does appear in NGC 4151, but only at wavelengths longer than about 5100 A. 
In the case of the blazar PKS 2155-304, the FR/RSS method confirms the ~2-hr X-ray/UV lag 
at a high level of confidence. 
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4. Discussion 

The results presented here establish the existence of wavelength-dependent continuum lags in 
at least one case. But do the measured lags indeed represent processes that are occurring in the 
continuum source itself? On account of the large line widths in the rich emission-line spectra of 
AGNs, there are essentially no truly "line- free" continuum windows; all wavelengths are at least 
somewhat contaminated by line emission. Is it possible, in fact, that these time delays are due to 
emission-line contamination? Excluding the contamination by the known strong emission lines, 
the wavelength-dependence of the lags in NGC 7469 (Fig. 6 of Wanders et al. 1997 and Figs. 
6 and 7 of Collier et al. 1998) argue fairly strongly against contamination by weaker emission 
features. Any possible emission-line contaminants would have to show a correlation of lag with 
wavelength that would reproduce the approximate r oc A 4 / 3 relationship apparently detected in 
NGC 7469 (Collier et al. 1998). It remains possible that the interband lags we have detected are 
due to unidentified contaminants, but we regard this as unlikely and will therefore assume that 
the measured lags are indeed a property of the continuum source. 

If the detected wavelength-dependent continuum lags arise as a consequence of the temperature 
structure of an accretion-disk continuum source, we need to ask whether or not the detection of 
wavelength-dependent continuum lags in NGC 7469 is consistent with the absence of detectable 
lags in NGC 3783, NGC 5548, and Fairall 9 and with the marginal detection in NGC 4151. Thus 
far we have concentrated only on the differences in the observations, and have not considered the 
differing intrinsic properties of these various galaxies and how this might influence the detection of 
such lags. In this section, we will assume that a simple thermal accretion disk model can account 
for the UV/optical continua of non-blazar AGNs, and we will scale the results for NGC 7469 to 
other AGNs to determine whether or not we should have expected to detect wavelength-dependent 
continuum lags in the programs that have been carried out. 

If we assume that the continuum is thermal emission from an accretion disk, and we 
further assume that the variability in different parts of the accretion disk is coupled by radiative 
reprocessing, then we find that the continuum at wavelength A is time-delayed relative to the 
primary source of photons such that 

r oc (MM) 1/3 A 4/3 , (1) 

where M is the AGN mass and M is the accretion rate (see Collier et al. 1998). We further assume 
that the UV luminosity is proportional to the accretion rate (i.e., we assume a constant efficiency 
for conversion of mass into light), and that the mass of the central source can be estimated virially 
from the broad-line widths ti FWHM and emission-line lags r LT , i.e., M oc Vp WUM T LT . By making these 
simple assumptions and scaling the results relative to NGC 7469, we can compute the expected 
UV-to-optical lag At under the assumption that M scales with the UV luminosity and that M 
scales like the square of the Lya line width times its time delay. We show the results of this naive 
calculation in Table 9. Column (1) gives the galaxy name, and columns (2) and (3) give M and 
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M, respectively, scaled to NGC 7469 assuming the relationships given above. These are used to 
compute a time delay At between the UV and optical wavelengths given in columns (4) and (5), 
respectively. This predicted time delay, scaled to the 1.8-day delay between the 1315 A and 6962 A 
variations in NGC 7469, is given in column (7). The expected lag given in column (7) should 
be compared with the 90% confidence interval (two-sided) on T cen t computed by the FR/RSS 
method, which is given in column (8). In the case of NGC 4151, the 90% confidence interval for 
the 1275 A/6925 A cross-correlation is surprisingly broad compared to the ±At68 width of only 
0.623 days (for the centroid). This is because of the highly non-normal CCPD, as shown in Fig. 4. 
In this case, the CCPD has a extended low-level tail, which is a characteristic of CCPDs based on 
a small number of data points (there are only 9 points in the 6925 A light curve). 

In each of the cases studied here, we find that the non-detection (or very marginal detection in 
the case of NGC 4151) of wavelength-dependent continuum lags is consistent with the predictions 
of the simple scaling model. The lags predicted by the scaling adopted are simply too small to 
have been detected by the experiments that have been carried out. Both the predicted lags based 
on scaling relative to NGC 7469 and zero lag lie safely within the 90% confidence interval in each 
case. 

Given this naive model and the considerable uncertainties in the determination of both M 
and M, we do not believe that previous non-detection of wavelength-dependent continuum lags 
poses any new theoretical challenge to our understanding of accretion disks in AGNs. 

5. Conclusions 

In this paper, we have introduced a model-independent Monte-Carlo method for estimating 
conservatively uncertainties in cross-correlation lags, and have tested this method for a variety of 
conditions. We have applied this method to the UV/optical light curves of several well-studied 
non-blazar AGNs and the blazar PKS 2155— 304. The reality of wavelength-dependent continuum 
lags found in NGC 7469 (Wanders et al. 1997; Collier et al. 1998) is supported by this analysis at 
at level of confidence higher than 95%. We confirm at about a 95% confidence level the detection of 
a '--'2-hr X-ray/UV lag in PKS 2155— 304. However, we find no evidence for statistically significant 
UV/optical continuum lags in the non-blazar AGNs NGC 5548, NGC 3783, and Fairall 9, and 
we find only marginal evidence for the effect in NGC 4151. By scaling other AGNs relative to 
NGC 7469, we find however that the absence of significant wavelength-dependent continuum lags 
in other AGNs does not appear to be inconsistent with a thermal accretion-disk model — most 
AGNs have not been monitored at high enough temporal frequency for UV/optical continuum lags 
to be detectable. Failure to detect such lags in previous experiments does not pose a serious threat 
to accretion-disk models, as long as the variability signal is propagated at the speed of light. 

We are grateful for support of this study by NASA (through ADP grant NAG5-3497 and 
LTSA grant NAG5-3233) and by the National Science Foundation (through grant AST94-20080). 
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Fig. 1. — Cross-correlation peak distributions for an FR/RSS simulation. The specific example 
shown is Test 2 in Table 2, i.e., the response of a thin spherical shell of radius R = 20 light days. 
Note that the distribution of CCFs peak values (r pea k, upper panel) is broader than the distribution 
for the CCF centroid (r cerit , lower panel). We find this to be generally the case. 



Fig. 2. — The UV time lags as a function of wavelength for NGC4151. The upper panel shows the 
average SWP (shorter wavelengths) spectrum and the average LWP (longer- wavelength spectrum). 
The second panel shows the root-mean-square (RMS) spectrum, in which the most variable parts 
of the spectrum stand out. The third panel shows the maximum correlation coefficient of the CCF 
fmax for each wavelength bin, and the bottom panel shows the CCF centroid T cen t of the CCF that 
results from cross-correlating the measurements at each particular wavelength with respect to the 
1300 A light curve. In the bottom panel, the highest peaks are due to the relatively slow response of 
LyaA1216 and CivA1549. Other than these features, there is no apparent wavelength-dependent 
trend as there was in NGC 7469 (Wanders et al. 1997). 



Fig. 3. — The optical time lags as a function of wavelength for NGC 4151. The upper panel shows 
the average blue and red spectra obtained at Wise Observatory. The second panel shows the 
RMS spectra computed from these same data. The light curves for each wavelength bin are cross- 
correlated with the 5125 A light curve; the third panel shows the peak value of the cross-correlation 
coefficient r max as a function of wavelength, and the centroid T cen t for each wavelength is shown 
in the bottom panel. The highest peak in the bottom panel is due to the slow response of the Ha 
emission line. There is also a hint of the He n A4686 response. The bottom panel shows that T cen t 
appears to increase with wavelength longward of about 5100 A. This appears to be similar to what 
was seen in NGC 7469. 



Fig. 4. — Cross-correlation peak distributions for the 1275 A/6925 A continuum bands in NGC 
4151. The CCPD shows a strong peak at at time delay of about 2.5 days, but slightly more than 
10% of the realizations yield values that form an extended tail towards lags smaller than about 2 
days. This occurs on account of the very small number (9) of data points in the 6925 A light curve; 
with so few points, it is generally not possible to determine a lag at a high level of confidence. 
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Table 1. A Typical Test 



Parameter 
(1) 


MC 
(2) 


FR/RSS 
(3) 


RSS 

(4) 


FR 

(5) 


T cent (days) 
Tpeak (days) 


19.11JI 
20l 2 


21.01J3 
20^ 


21-0iy 
20^ 


21.0113 
201^ 



Table 2. Repeated Tests with Constant Parameters — Thin-Shell Model 



Test 


Tcent 


(days) 


Tpeak (days) 


No. 


MC 


FR/RSS 


MC 


FR/RSS 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


i9.ni:| 


21.0113 


2011 


2011 


2 


i7.i+l;l 


16 2+ 4 - 5 
i0 - z -2.6 


2011 


24+ 1 


3 


i9.5+l:I 


23.9+^ 


2011 


24+ 2 


4 


20.0lj| 


18.411° 


20+1 


18+1 


5 


18.0113 


18.7+J3 


211J 


2111 


6 


20.2+H 


24.01U 


20+1 


22+1 


7 


18.9113 


16.HII 


19+| 


16+? 


8 


19.1111 


19.6+i| 


19+1 


21+ 2 3 


9 


99 Q+3-9 


17.611"? 


2111 


19 + 5 


10 


18.5l^ 7 5 


20-0111 


201^ 


2111 
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Table 3. Repeated Tests with Constant Parameters — Thick-Shell Model 



lest 


Tcent (days) 


Tpeafc (days) 


AT „ 

1NO. 


MC 


FR/RSS 


MC 


FR/RSS 


(i) 


(2) 


(3) 


(4) 


(5) 


i 


ii-56l8:S 


10.7711:11 


i2.oil:8 


n.oil:8 


2 


10.3ll|;| 


9.55_ 1 24 


n.o+18 


io.sll-8 


3 


13.85_ L66 


11.76l?-g 


12.0 + f8 


14 +0 - 5 


4 


10.801°$ 


ii.ool? : | 


n.oil:8 


12.018:8 


5 


12-1911:11 


12.38_ 3 03 


i2.o!£{j 


i2.5l 3 :8 


6 


9-96±i:90 


9.87l|| 


ii.ol|8 


11-5118 


7 


8.98+1-lf 


9.75l :l 


io.ol£|j 


9.0+fg 


8 


9.56ll:| 


8.5811:1 


ii.oiI:8 


9.ol? : g 


9 


6.42ll- 8 7 ° 


9.051^1 


8-5118 


8.5i?:8 


10 


11.8518U 


19 oc;+ 1 - 67 
iz.zd_ g4 


ii.5ll:8 


ii.ol8:i 



Table 4. Thin-Shell Models in the Poor-Sampling Limit 



PDS R T cent (days) r peak (days) 

index (It days) N MC FR/RSS MC FR/RSS 
(1) (2) (3) (4) (5) (6) (7) 



-2.5 


20 


40 


i8.5ii8 


14.818:1 


1911 


1911 


-2.5 


20 


20 


i9.8±f:f 


17.8l8i 


19+1 


1918 


-2.5 


20 


10 


19.9ll| 


i8.6i 6 7 :8 


1911 


2018 


-2.5 


2 


40 


2 +a9 
z - u -i.o 


i.9l?:8 


211 


2 +1 


-2.5 


2 


20 


i.9ll:8 


1.51^ 


211 


2 +1 
z -3 


-2.5 


2 


10 


o.sllJ 


°-'-15.7 


lit 


4 -20 


-1.8 


20 


40 


i8.3±| : g 




19+? 


22 + s 


-1.8 


20 


20 


19.4111 


18.0111 


2118 


1411 1 



-1.8 20 10 18.5lf;| l 7 - 6 ^ 7 19l? 20l^ 7 
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Table 5. NGC 7469 Cross-Correlation Results 



Centroid Peak 



First 


Second 


Tcent 




Tpeak 




Series 


Series 


(days) 


P(< Odays) 


(days) 


P(< Odays) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1315 A 
1315 A 
1315 A 


1825 A 
4845 A 
6962 A 


o.22±g : g 

-i oc+0.48 
i ' zo -0.35 
-, O/I+0.93 
i -° 4 -0.94 


0.024 
0.004 
0.030 


0.073;" 


0.166 
0.006 
0.046 



Table 6. UV/Optical Continuum Lags for Well-Studied AGNs 

Centroid Peak 

Galaxy (days) P(< Odays) (days) P(< Odays) 

(1) (2) (3) ' (4) (5) ' 



NGC 5548 WE 1989 a 


1 3+ L5 


0.351 


9 r+1.5 
z - J -3.0 


0.205 


NGC 5548 HST/IUE 1993 b 


i O+0.6 
L - z -1.2 


0.191 


o.7±i:| 


0.341 


NGC 3783 c 


0.2!! 2 , 


0.478 


0.9+?; 8 7 


0.354 


Fairall 9 d 


-2-8+i;| 


0.622 


-0-510 


0.584 



a UV: 1337 A flux (GEX version) from Clavel et al. (1991), Table 2. Optical: 
5100 A flux from Peterson et al. (1992), Table 12. 

b From Korista et al. (1995). UV: 1350 A flux from Table 12, plus selected IUE 
SWP fluxes before JD 2449079.0 from Table 16. Optical: 5100 A fluxes from 
Table 21. 

C UV: 1460 A flux (GEX version) from Reichert et al. (1994), Table 5. Optical: 
5150 A flux from Stirpe et al. (1994), Table 7. 

d UV: 1390 A flux from Rodrfguez-Pascual et al. (1997), Table 2. Optical: 
5340 A flux from Santos-Lleo et al. (1997), Table 2. 
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Tablc 7. NGC 4151 Cross-Correlation Results 









Peak 




Centroid 


First 


Second 


Tpeak ± AtmL & 


Tpeak 




Tcent 




Series 


Series 


(days) 


(days) 


P(< Odays) 


(days) 


P(< Odays) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1275 A 
1275 A 
1275 A 
1275 A 


1820 A 
2688 A 
5125 A 
6925 A 


-o-oi±8:8| 
o.ooi8:g? 
-o.o5^:| 


-0.02518:83 

o ooo+ - 095 

u.uuu_ 00 5 

U.4WU_ 1 2 6 9 

9 /I8^+ - 280 
/.48t>_ 0450 


0.930 
0.470 
0.449 
0.129 


-o.o49i8:8l 
o.o86i8:8I 

0.2961S 
2.46218;^ 


0.991 
0.042 
0.423 
0.102 



a From Edelson et al. (1996) Table 4. 



Table 8. PKS 2155-304 Cross-Correlation Results 



Centroid Peak 



First 
Series 
(1) 


Second 
Series 
(2) 


Tcent 

(hrs) 
(3) 


P(< Ohrs) 
(4) 


Tpeak 

(hrs) 
(5) 


P(< Ohrs) 
(6) 


24 A 
1400 A 
2800 A 


1400 A 
2800 A 
5000 A 


o.o5ll:l 
-o.7ol 2 :^ 2 


0.052 
0.520 
0.729 


o 9+1.0 
z - z -1.6 

o.iil: 2 
-o.6li3 


0.034 
0.621 
0.676 



a Confidence levels are 90% for direct comparison with Edelson et al. 
(1995) Table 2. 
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Table 9. Comparison with Theoretical Scaling 









Ai 


A 2 


At 


90% Confidence 


Galaxy 


M a 


M a 


(A) 


(A) 


(days) 


Interval (days) 


(1) 


(2) 


(3) 


(4) 


(5) 


(7) 


(8) 


NGC 4151 


0.53 


0.44 


1275 


6925 


1.1 


[-1.90, +3.00] 








1275 


5125 


0.7 


[-1.25, +1.40] 


NGC 3783 


2.92 


0.38 


1460 


5150 


1.1 


[-4.69, +4.62] 


NGC 5548 


2.00 


1.04 


1350 


5150 


1.4 


[-0.54, +2.49] 


Fairall 9 


4.74 


8.42 


1390 


5340 


4.0 


[-9.63, +10.28] 



a Relative to NGC 7469. 
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